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We propose a new model in order to study behaviors of self-organized system such as a group 
of animals. We assume that the individuals have two degrees of freedom corresponding one to 
their internal state and the other to their external state. The external state is characterized by 
its moving orientation. The rule of the interaction between the individuals is determined by the 
internal state which can be either in the non-excited state or in the excited state. The system is put 
under a source of external perturbation called “noise”. To study the behavior of the model with 
varying noise, we use the Monte-Carlo simulation technique. The result clearly shows two first-order 
transitions separating the system into three phases: with increasing noise, the system undergoes 
a phase transition from a dilute disordered phase to an ordered compact phase and then to the 
disordered dispersed phase. These phases correspond to behaviors of animals: uncollected state at 
low noise, flocking at medium noise and runaway at high noise, respectively. 


I. INTRODUCTION 

The collective behavior of animals is a widely observed 
phenomenon in various biological systems. It is one of 
the main topics which have been extensively investigated 
during the last two decades using methods borrowed from 
many different areas of science including physics, ap¬ 
plied mathematics and engineering. For recent papers, 
the reader is referred to works of Couzin and cowork¬ 
ers, of Albano and his group among others [m- In 
these works, many aspects of animal groups have been 
discussed, among which the collective motion and their 
origin, the relation between individual and collective be¬ 
haviors and the relation between group size and collective 
decision-making. 

The flocking is a behavior of some animal species where 
they stay together in a group for social reasons. They de¬ 
rive many benefits from this behavior including defence 
against predators, easier collective moving, enhanced for¬ 
aging success and higher success in finding a mate. When 
they are faced with a danger such as predators, their nat¬ 
ural instinct is to flee not to fight. They use their natural 
herding instinct to bind together in a group for safety. All 
individuals of the group will move away from the preda¬ 
tor in the same direction and then stampede as fast as 
they can when being under the predator’s attack. If there 
is no danger, then they are spread to find foods instead 
of staying in the flocking state. Many experimental facts 
and observations have been for example mentioned in the 
review of Vicsek and Zafeiris [I| . Well-known examples 
are found in populations such as large schools of fish Q 
or gatherings of birds Q . Biologically, it is known that 
the flocking behavior is advantageous for survival of a 
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population [^-[TT| : reducing the risk of capture by preda¬ 
tors, increasing higher mating efficiency, easier search for 
food, efficient learnin g of external stimuli, and reducing 
overall aggression |12 h1 ^ . 

In 1987, Reynolds first suggested a simple model con¬ 
sisting of three rules: separation, alignment, and cohe¬ 
sion rules [l^ . These rules describe the behavior of each 
individual in interaction with other neighboring individ¬ 
uals. All or some of the three rules were mathematically 
expressed and then analyzed by Vicsek and his cowork¬ 
ers They mainly focused on the transition 

between coherently moving and runaway in a stampede. 

The flocking behavior has been conventionally stud¬ 
ied through simulation in two frameworks: population 
(Eulerian or continuum models) and individuals (agents 
or particle-based models) [a la [2l|. In the popula¬ 
tion framework, the flock was collectively addressed while 
flock-density was used as a key variable to present spatial 
and temporal dynamics of aggregation frequently with 
partial differential equations of advection-diffusion reac¬ 
tion [ 23 . In the individual framework, the flock of 

agents has been simulated by using ordinary and stochas¬ 
tic equations of motion to describe interactions among 
agents 0, [13) [23 - 1^ . This approach attempted to repli¬ 
cate naturally observed phenomena from not only animal 
groups but also other self-propelled characteristics [l^l 
and to compare the evolved characteristics with those of 
actual animal flocking [2^ in order to better understand 
possible mechanisms by which these characteristics may 
have evolved. 

The previous models have predicted with success flock¬ 
ing behavior of animal groups at high noise. They have, 
however, difficulties to explain uncollected states at very 
low noise. This has motivated our present work: in addi¬ 
tion to the interaction between neighboring animals, we 
introduce an internal degree of freedom which indicates 
whether an animal is in a non-excited or an excited state. 
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we will show that the collective behavior of animals in the 
whole range of noise can be explained. 

Section II is devoted to the description of the model. 
Section III shows the Monte Carlo (MC) simulation re¬ 
sults for testing our model. Conclusions are given in Sec. 
IV. 


II. THE MODEL 

In biology, all the members of an animal group are 
spread to find foods if there is no danger. In this situa¬ 
tion, they are distributed in the space with a small con¬ 
centration and out of alignment. So, we say the group 
of animals is in a “uncollected” behavior. When the ani¬ 
mals are faced with danger such as predators, they bind 
together in a small area for safety with the same orienta¬ 
tion and high concentration. This state is called “flock¬ 
ing” state. Facing a danger, animals will move away from 
the predator in the same direction and then stampede as 
fast as they can when being under the predator’s attack. 
At the hnal stage, they are in a “runaway” state. 

In order to study the phase transition behavior, we 
have to map the group of animals into a physical system. 
We consider an animal as a particle i with two degrees 
of freedom: an external parameter ai characterizing the 
animal orientation which depends on its interaction with 
the others, and an internal parameter Si indicating either 
it is in the non-excited {Si = 0) or in the excited {Si = I) 
individual state. 

Let us define the internal degree of freedom. The inter¬ 
nal state of a particle can be described by two levels: the 
first level is characterized by a negative energy — e which 
describes the calm, unworried stable state, and the sec¬ 
ond level by a positive energy -\-e which expresses some 
degree of anxiety. Introducing an external noise 77 which 
plays the role of the temperature in statistical physics, we 
can write the canonical probability of these two states as 


e/77 


3-e/’; 


Pi = 


P2 = 


z z 

where 0 is the partition function given by 
2 _ gV'/ _|_ — 2 cosh(e/? 7 ) 


( 1 ) 


( 2 ) 


Note that pi > p 2 because e > 0. So, the “positive” ten¬ 
dency towards tranquility is proportional to pi—p 2 - The 
total number of particles having this tranquility state 
which is called non-excited state is thus given by 


No = N{pi -P 2 ) = N 
= iVtanh(e/ 77 ) 




( 3 ) 


where N is the total number of particles in the system. 
The number of excited particles is then 


N, = N - N, 


0 


( 4 ) 


One sees in the above equations that in zero noise one has 
Nq = N, namely there is as expected no excited particles. 
While, when the noise goes to infinity, one has Nq = 0 
and Ne = N: all particles are excited. At a given noise, 
statistically one has a well-defined Nq and Ne with the 
conservation N = Nq + Ne. 

In the simulation which will be done below at a given 
noise, we generate at random Nq according to Eq. @. 
The remaining particles in the system are thus in the 
excited state. 

Now, we put this system of individuals on the lattice 
where each individual can move on 2D triangular lattice 
of linear size L. The number of lattice sites should be 
greater than that of individuals, i.e., ^ N. We denote 

by tJi the orientation of individual Sp. Ui is dehned as in a 
g-state Potts model, i.e. ai = 1, 2,..., g. For simplicity, 
we consider g = 6, so the orientations <77 = 1 , 2 ,..., 6 of 
individual Si can be defined by the vectors which connect 
a site to its nearest neighbors (NN) with the following 
angles measured from the x axis: ipi = 0, 7r/3,..., Stt/S. 
The interaction between nearest-neighboring animals is 
given by the Hamiltonian 

H = ^ cos[Tr{a^ - aj)/3], (5) 

<i,j> 

where the sum J2<i j> made over up to the third near¬ 
est neighboring individuals Si and Sj. For simplicity, 
Kij is assumed to take the form of the Lennard-Jones 
potential: 


= AJ [{ro/n^jf^ - (ro/r,j)®] , 

where Vi^ is the distance between two individuals. We 
choose ro = 0.89 in order that Ki^j ~ J at Vij = 1 which 
is the lattice spacing, J being a positive constant. The 
interaction between two individuals depends on their in¬ 
ternal state: Kij = 0 if both individuals are non-excited 
Si = Sj = 0, and Kij y/ 0 if otherwise. We empha¬ 
size that the interaction rule we impose means that an 
excited particle can interact with non-excited particles, 
dragging them into a collective motion. 

Let us note that, without the internal degree of free¬ 
dom and without particle motions, the above Hamilto¬ 
nian is the localized Potts clock model which has been 
solved for g = 2, 3 and 4 [s^ . The case of very large g has 
been solved by Frohlich and Spencer [3l|. However, for 
g = 5, 6 , ... there are not (yet) exact results. Of course, 
our model is more complicated because the particles are 
mobile with an internal degree of freedom. 


III. THE PHASE-TRANSITION BEHAVIOR 

We use the MC simulation technique to study the 
above model. The main physical quantities such as the 
order parameter Q and the concentration p are dehned 
in what follows. 

For Q, we have 
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where 


with 


Q 


g-^max ~ 1 


-Mmax 


max(Mi, M 2 , Mq) 

N 


Mi — {i — 1,2,q) 
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( 6 ) 

(7) 

( 8 ) 


where the sum on j is performed over all sites of the 
system. One sees that in the ordered state where there 
is only one kind of aj one has Q = 1. In the disordered 
phase where all values of aj are equally present, namely 
Mi = N/q for any i = 1,..., q, one has Q = 0. 

For the concentration p we have 


1 . ' 

( 9 ) 

i=l 

where is the number of NN individuals around Si. The 
quantity p in Eq. m characterizes the spatial distribu¬ 
tion of the individuals. The behavior of animals can be 
adequately described by the two parameters Q and p. 

In the simulations, we use e = 0.04, J = 1.0 (taken as 
the unit of energy) and N = 100, 400 and 900 with the 
lattice size L = 4 x ^/N. For the initial configuration of 
the system, we randomly generate the particle positions 
on the lattice and take orientations of all individuals from 
an uniform distribution. At each MC step, we randomly 
choose Nq individuals (according to Eq. ([3])), and set 
them to be non-excited Si = Q. The remaining particles 
are in the excited state Si = 1. Their position and ori¬ 
entation are updated by Metropolis algorithm. At each 
rj, the equilibration time lies around 4 x 10® MC steps 
per individual and we compute statistical averages over 
8 X 10® MC steps per individual. Each time a particle 
gets out of the lattice at one boundary, we put it into 
the system by the other end using the periodic bound¬ 
ary conditions in order to conserve the total number of 
particles. 

We plot in figures [1] and [3] the order parameter and 
the concentration as a function of external noise for sev¬ 
eral system sizes. It clearly shows the existence of three 
phases which are separated by the two transitions at very 
low noise and high noise. In phase I at low noise, the 
system is in the disordered phase with low concentra¬ 
tion (Q Ri 0 and p r: 0). This phase is equivalent to 
the uncollected behavior of an animal group which is dis¬ 
persed over the whole space because almost of them are 
non-excited, it is called “free” phase with a few “con¬ 
tacts” between the individuals. Phase III, at high noise, 
is spatially sparse as the phase I, but the particles are 
very mobile in all directions. It describes the runaway 



FIG. 1: Order parameter versus noise with the system sizes 
N — 100 (circles), 400 (squares) and 900 (diamonds). The 
insets show the enlarged scale at low (a) and high (b) noise. 



FIG. 2: Goncentration versus noise with N = 100 (circles), 
400 (squares) and 900 (diamonds). 


behavior. At medium noise, the animals are in phase II 
where they are compactly moving in an ordered phase 
with Q R 1 and the maximum concentration p r 6. This 
phase corresponds to the flocking behavior. 

The spatial distribution of the individuals is also pre¬ 
sented in Fig. [3] with snapshots taken at several values of 
noise rj = 0.0522,0.1367,0.5000 and 0.5667, at the final 
stage of a MC simulation. The vectors indicate their po¬ 
sition and moving orientation. One sees that Figs[3Ka,d) 
show that the system has the same distribution where 
they are randomly distributed in the plane with differ¬ 
ent orientations. However, the instantaneous snapshots 
do not permit to see the difference between two phases: 
the particles in phase I are almost immobile while they 
are moving very fast in a disordered manner in phase III 
(this difference is seen in videos). Fig. [Sjb) shows the 
flocking behavior, and Fig. [SKc) shows the distribution 
of the system at the noise closes to the II-III transition 
point. 

At this stage, let us examine again Fig. [T] it shows 
discontinuities of the order parameter at both transitions 
I-II and II-III, indicating a signature of a first-order tran¬ 
sition. This abrupt behavior corresponds to the fact that 
the animals immediately flock when they are faced with 
danger, and they runaway when being under a predator’s 
attack. These behaviors are more likely a first-order tran- 
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FIG. 3: Snapshots at different values of noise with N = 100. 
The arrow indicates the orientation of individual. See text for 
comments. 



FIG. 4: Histogram versus energy at low (a) and high (b) 
noise at the transition points, with the system sizes N — 100 
(solid), 400 (dotted) and 900 (dashed). 


sition. To confirm the first-order transition we have per¬ 
formed the energy histogram at the transition point. We 
show in Fig. 0] the energy histogram at two transitions 
I-II and II-III, for three system sizes. The double-peak 
histograms are clearly shown for both transitions at low 
noise (a) and high noise (b). The dip between the two 
maxima becomes deeper with increasing size. Note that 



FIG. 5: Order parameter versus noise with N — 100 for the 
localized lattice model (diamonds) and lattice gas model (cir¬ 
cles). 


a “true” discontinuity happens only when the dip comes 
down to zero. The distance between the two peaks is 
then the latent heat. To see this, we need sizes much 
larger than N = 900. But this is out of the scope of our 
present purpose. 

Let us discuss on the phase transition at high noise. 
It is known that the g-state Potts model localized on 
a lattice has a first-order transition with g > 4 in two 
dimensions. The clock Potts model localized on a lat¬ 
tice was solved oi^ for g < 4 [3^, though it was solved 
for very large g [3l|. The model treated in the present 
paper corresponds to a mobile Potts model where par¬ 
ticles can go from one lattice site to another. This was 
never solved before. Note however that the transition II- 
III corresponds not only to an orientational disordering 
of the Potts parameter but also to the breaking of the 
compactness of the flocking state. This is similar to a 
melting transition where the solid phase melts into the 
liquid phase. The melting has often a first-order char¬ 
acter in three dimensions [32|. In two dimensions, it is 
known that long-range solid ordering does not survive 
at finite temperatures if the system has isotropic short- 
range interactions according to Nelson and Halperin [^ : 
at a first critical temperature, bound pairs of dislocations 
formed at low temperatures are unbound, giving rise to 
a phase with no translational ordering but with orien¬ 
tational hexatic structure. The latter phase undergoes a 
“critical” phase transition to the disordered phase. It was 
however found by MC simulation that the second transi¬ 
tion is not critical but it is a first-order two-dimensional 
disclination melting 3^. Our results shown in Figs. [3] 
(c,d) and 0] indicate that our model undergoes a first- 
order transition from the hexatic ordered phase to the 
liquid phase, in agreement with earlier MC simulations 


Finally, to compare our mobile Potts model with the 
corresponding localized Potts model, we have simulated 
the same system defined by m including the internal 
degree of freedom on a lattice of size = N, i.e. the 
individuals are not allowed to move on the lattice. We 
show in Fig. 0] results of both mobile and localized Potts 
models: the transition temperature is about 0.5844 for 
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the mobile case, while it is ry ~ 1.1126 for the localized 
Potts model. 


IV. CONCLUSIONS 

We proposed in this paper a simple model for studying 
the behavior of animal groups as a function of an exter¬ 
nal perturbation, called “noise”, such as dangers coming 
from predators. We showed that with increasing noise, 
the system has three phases I, II and III separated by 
two transitions, the first transition occurs at a low noise 
and and the second one at high noise. The three phases 
are disordered, ordered and melted phases which corre¬ 
spond respectively to the following behaviors of animals: 


uncollected state, flocking state and runaway state. Both 
transitions from one phase to another are found to be of 
the first-order. 

This model is similar to the lattice gas model devel¬ 
oped by Csahok et. al. [l^ but in our model the align¬ 
ment rules of the individuals depend on the internal state, 
excited or non-excited, of its neighbors including itself. 
Therefore, our model can be applied not only for study¬ 
ing the phase transition behavior of animals at high noise 
as in previous models, but also for analyzing the animal 
behavior at low noise where animals are in a uncollected, 
disordered phase. 

This work was supported by the Nafosted (Vietnam 
National Foundation for Science and Technology Devel¬ 
opment), Grant No. 103.02-2011.55. 
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